ALBERTA THY-3-92 



ALL-ORDERS WORMHOLE VERTEX OPERATORS 
FROM THE WHEELER-DEWITT EQUATION 

ALEX LYONS* 

Theoretical Physics Institute, Department of Physics 
University of Alberta 

Edmonton, Alberta, T6G 2J1, CANADA 

Submitted to Nucl. Phys. B, Jan 1992 
Abstract 

We discuss the calculation of semi-classical wormhole vertex operators from 
wave functions which satisfy the Wheeler-deWitt equation and momentum con- 
straints, together with certain 'wormhole boundary conditions'. We consider a 
massless minimally coupled scalar field, initially in the spherically symmetric 
'mini-super space' approximation, and then in the 'midi-superspace' approxi- 
mation, where non-spherically symmetric perturbations are linearized about 
a spherically symmetric mini-superspace background. Our approach suggests 
that there are higher derivative corrections to the vertex operator from the 
non-spherically symmetric perturbations. This is compared directly with the 
approach based on complete wormhole solutions to the equations of motion 
where it has been claimed that the semi-classical vertex operator is exactly 
given by the lowest order term, to all orders in the size of the wormhole throat. 
Our results are also compared with the conformally coupled case. 
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1 Introduction 



Wormholes can be thought of as euclidean solutions to the field equations for gravity, 
possibly coupled to some matter system, which connect two asymptotically flat regions.* 
These solutions can be thought of as contributing to the semi-classical (zero-loop) approxi- 
mation to the partition function for gravity, and thereby to the Green functions for fields in 
the asymptotic regions. There are, however, problems with just considering real euclidean 
solutions. Firstly, it is known that a complex contour is needed in order to make the path 
integral for gravity well defined, and it is not clear how that contour should be defined, 
particularly when matter fields are present. Thus it would seem that arbitrary complex 
saddle-points, in which both geometry and matter are complex, might be just as relevant to 
the semi-classical evaluation of the partition function as real euclidean solutions. Secondly, 
it is known that in order for real euclidean wormhole solutions to exist, one needs very 
particular types of matter fields, such as a conformally coupled, or imaginary minimally 
coupled, scalar field. However, for wormholes to provide a viable solution to the cosmo- 
logical constant problem [2] , or to be relevant to the late stages of black hole evaporation 
[3], one would like their existence to be independent of the particular matter fields in the 
system. 

It is for these reasons that Hawking and Page [4] introduced the idea that wormholes 
should be considered instead as solutions to the Wheeler-deWitt equation. The wave func- 
tions of wormholes would satisfy particular boundary conditions corresponding to the fact 
that they describe states which, semi-classically, correspond to four-geometries connecting 
two asymptotically flat regions. Hawking and Page show that one can find a complete spec- 
trum of such wave functions which generalizes to a rather large class of matter field models. 
However, they did not extract a low-energy effective action from their wave functions. It is 
generally expected that at low energies it is possible to replace the wormhole ends by ver- 
tex operators in the effective field theory, which would then give rise to an effective bilocal 
action. The aim of this paper is to show that it is possible to extract the vertex operators 
directly from the wave functions. The approach based on wormhole wave functions turns 
out to be more general than using euclidean solutions to the field equations, in that it 
potentially allows a wider class of operators in the effective action. In general we would 
expect the effective action to include all possible gauge-invariant and Lorentz-covariant 
operators. The euclidean solution approach fails to find all these operators, due to the 
restricted nature of possible solutions. For example, in the imaginary scalar field model, 

By 'asymptotically flat' we mean 'asymptotically euclidean in the sense of Gibbons and Pope [1]'. 
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the euclidean wormhole solutions have to be O (4) -symmetric. Grinstein, Maharana and 
Sudarsky have calculated the semi-classical vertex operator corresponding to these worm- 
holes, to all orders in the wormhole throat radius [5], and they show that there are no 
derivative terms in the effective action. However, other work based on the wormhole wave 
functions suggests that one obtains these derivative terms [3, 6-8]. 

We would like to point out that considering just solutions to the Wheeler-deWitt 
equation is, itself, probably an approximation. It is now realized (using a simple two- 
dimensional model) [9] that the wormhole wave functions need not actually be 'on-shell', 
i.e. satisfying the Wheeler-deWitt equation. In the two-dimensional model considered in 
[9], the 'on-shell' states only provide the dominant contribution in the situation where 
the wormholes are long and thin. The 'on-shell' states appear as poles in the vertex 
operator, when the integration over the size of the wormhole throat is performed. In 
four-dimensional gravity we do not know the measure in the integration over the size of 
the wormhole. However, we shall assume that it is such that 'on-shell' states provide the 
dominant contribution to the low-energy effective action, as in the two-dimensional model. 

The plan of this paper is as follows: In Sec. 2, we consider the spatially homogeneous, 
O (4) -symmetric, massless scalar field model. We show how the : e lfc< ^: effective interactions, 
usually associated with the euclidean axionic (or imaginary scalar field) wormhole solutions, 
can be thought of as arising from a basis of solutions to the Wheeler-deWitt equation. 
Thus the wave function approach (in the O (4) -symmetric mini-superspace truncation) 
reproduces the results usually obtained by considering euclidean wormhole solutions. In 
Sec. 3, we discuss the non-spherically symmetric perturbations of the axion wormhole in 
the wave function approach. We treat the problem in the 'midi-super space approximation' 
as used in [10]. Excitations of the lowest spatially inhomogeneous mode (n = 2) would be 
expected to produce a tower of effective interactions of the form : (<9(/>) m e lfc< ^: , for even m. 
This is based on the fall-off behaviour of the modes in the asymptotic region. However, a 
more careful analysis of the linearized lapse and shift constraints shows that only the ground 
state wave function is allowed. This leads to just the vertex operator with m = 0. In 
contrast, for the higher spatially inhomogeneous scalar modes one obtains higher derivative 
corrections to the vertex operator. Our findings for the lowest mode are in agreement with 
the 'non-renormalization theorem' of Grinstein et al [5]. However, it is claimed in [5] that 
the vertex operator :e lk ^: is correct to all orders in the size of the wormhole. In our 
formalism, one would expect higher derivative corrections to this operator. We discuss 
the differences between our conclusions and those of [5] and also compare our results with 
those of Dowker [11] in the conformally coupled case. We summarize our main results in 
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Sec. 4. Finally, in the Appendix, we discuss in more detail the construction of the solutions 
to the euclidean Schrodinger equation which are used in Sec. 3. 

2 Wormholes, scalar fields, and the Wheeler-deWitt equation 

The aim of this section is to relate the wave function description of wormholes to 
the more usual description in terms of real euclidean wormhole solutions. We start by 
considering the simplest mini-superspace model which exhibits wormhole solutions. This is 
the model consisting of gravity minimally coupled to a massless scalar field. The scalar field 
can be thought of as the Goldstone boson corresponding to a spontaneously broken global 
U(l) symmetry [12]. Alternatively there is a dual description of this model (related by a 
simultaneous canonical and Hodge duality transformation) in terms of an anti-symmetric 
three-index field strength. It was in this form that wormhole solutions were first found 
[13]. However, we shall here (for simplicity) work with the scalar field formulation. In this 
model the spatial sections are taken to be three-spheres and the metric ansatz takes the 
form: 

ds 2 = a 2 (N 2 (t)dt 2 + a 2 (t)dn 2 ) , (2.1) 

with the matter field taken to be homogeneous on the spatial sections. The constant 
a 2 = 2/3nm 2 is included for later convenience, and dfl 2 is the usual round metric on the 
unit three- sphere. 

With this ansatz the euclidean action for this mini-superspace model is 

This is the action appropriate for the boundary value problem where a and <p are fixed on 
the two boundaries, labelled by t± . The boundary term in the action is present in the case 
where the boundary labelled by t- is taken to become asymptotically flat space. This is 
the case which we shall be interested in throughout this paper. The limit t- — > — oo is the 
asymptotically flat region of the space-time. Varying the action with respect to N and <p 
leads to the constraint equation (the hamiltonian constraint for gravity) and the equation 
of motion for which are: 

PPVdt) ~ + N 2 \dt) 

V 7 V 7 (2.3) 

d_ (a?df\ 
dt \N dt J 
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These form a complete set of equations for the system, since if these hold at all times then 
the a equation of motion is automatically satisfied. 

The simplest general form of the solution to these equations can be obtained in the 
gauge where N = 1/a. In this case (2.3) can be integrated trivially to yield the solution 

ds 2 = a 2 (a- 2 (t)dt 2 + a 2 (t)dQ 2 3 ) 
a 2 = v 7 ^ 2 - Q 2 



1, (2t-Q 
- m 



(2.4) 

where the constant of the motion (the euclidean momentum conjugate to 4>) is 

Q = a 4# =a 2 sinh2 (0-0_) . (2.5) 

Here </>_ is the asymptotic value of 4> at infinity, and a constant of integration can also be 
added to t . 

Flat space is recovered as the special case (5 = 0, while the wormhole solution is found 
by taking Q to be purely imaginary. The metric then has two asymptotically flat regions, 
as t — > ±oo , with a throat at t = . This solution will also have a purely imaginary . 
One might argue that this solution is a bit of a cheat, since the matter field has negative 
definite energy density and seems rather unphysical. However, there are two good reasons 
for arguing that such a configuration may have relevance. Firstly, from the point of view 
of providing a semi-classical contribution to the partition function, it is known that the 
contour of integration must be distorted into the complex plane for convergence of the 
path integral. In that case one might argue that any saddle-point may be important, even 
one in which the metric and matter field become complex [14]. The second argument 
comes from viewing the semi- wormhole solution (— oo < t < 0) as providing a semi- 
classical contribution to the amplitude for a tunnelling process in quantum gravity. The 
idea is that if one wants the amplitude to tunnel between states of definite real lorentzian 
momentum for the scalar field, then the boundary data ensures that Q must be purely 
imaginary, as in this solution [15]. 

In the alternative treatment of wormholes which uses the Wheeler-deWitt equation, 
we can bypass the argument over whether or not one should consider complex solutions. 
One might expect that a WKB approximation to a solution of the Wheeler-deWitt equa- 
tion could be written in the form e ~ Isp ^ a '^ , where I sp (a,4>) is the euclidean action of 
the classical solution (2.4) between asymptotically flat space as f ^ — oo, and an inner 
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boundary given by t + , on which a and <fi are fixed and real. The action of this classical 
solution, given by substituting (2.4) into (2.2) , is 



I sp {a,4>) = -t+ = - 



a 2 cosh 2(0 - 4>_) 



(2.6) 



which reproduces the flat space action in the case </> = </>_. Notice that the solution which 
we use has real <f>. We do not need to go to complex 4>, because we do not need to find a 
classical solution which has two asymptotically flat regions. Instead, we are interested in 
the wave function for real values of its arguments. The Wheeler-deWitt equation for this 
model is 



where the factor ordering has been chosen to be the one which is covariant under general 
co-ordinate transformations in the (a, 4>) mini-superspace. The family of wave functions 



in fact solve this equation exactly, although the derivation from the euclidean action of a 
classical solution suggests that one would have only expected to satisfy this equation to 
leading order in the WKB approximation. 

The following questions now arise: Firstly, in what way are the wave functions i/V- 
related to the euclidean wormhole solution to the field equations with an imaginary scalar 
field? If there is a straightforward connection between them, one can then ask whether 
the low-energy vertex operators for these wormholes can be obtained directly from the 
wave functions. In the past vertex operators have been obtained by considering Green 
functions on the euclidean background solution [16]. One motivation for using the wave 
functions instead comes from the possibility of generalizing this model to include spatially 
inhomogeneous perturbations. One expects that these extra degrees of freedom would 
enlarge the space of wormhole quantum states. The vertex operators corresponding to 
these extra states would perhaps be missed in a treatment which just considers the field 
theory on a background 0(4) symmetric solution. These will be the questions we shall 
consider in this paper. 

To address the first question, we first consider the Fourier transform (with respect to 
4>) of our wave function. In other words we consider 




(2.7) 



±a 2 cosh2(<£-<£_) 



(2.8) 




(2.9) 



5 



This can be thought of as changing from the (a, 0) to the (a, k) representation of the state 
. t(j is given explicitly by a modified Bessel function of imaginary order, Ki ik (a 2 /2) , 
multiplied by e lk< ^- . It has the approximately exponential form [17] 



^~ (a 4 -£; 2 )- 1/4 exp 



— - a 4 — k 2 + k arcsin(/c/a 2 ) j + ik(p- 



(2-10) 



in the region a 4 3> /c 2 and /c 2 3> 1 . 

This form is suggestive of a euclidean saddle-point approximation to a path inte- 
gral representation for if) . Indeed, consider the complex canonical transformation from 
(a(t), <j>(t)) to (a(£), k(t)) , where fe = — ia 3 0/iV. For a 4 > /c 2 , the classical variational 
problem with (a, fc) fixed on the inner boundary and — > 0_ at infinity has a euclidean 
stationary point whose action is precisely the negative of the value in the exponent above. 
The action (2.2) has to be supplemented by a matter field dependent boundary term on 
the inner boundary given by — ik(t+)(j)(t+) , since the classical variational problem has 
changed from one in which <p is fixed on the boundary to one in which k is fixed. This 
boundary term is fixed by the mathematical requirement that the variational principle 
subject to the new boundary conditions should yield the equations of motion. It is this 
boundary term which yields the \k arcsin(/c/a 2 ) — ikcj)- term in the evaluated euclidean 
action. The \\/ a 4 — k 2 term comes from (2.2) evaluated at the stationary point. 

The saddle-point four-geometry in the (a, k) boundary value problem is precisely a 
part of the imaginary scalar field wormhole outside a three-sphere of given radius aa , car- 
rying 'charge' k. (This would be the lorentzian U(l) charge carried by the wormhole in 
the Goldstone boson interpretation of the theory.) The boundary value problem has two 
solutions however, unlike in the case where the field <p is specified on the boundary. These 
correspond to taking more than or less than half of the wormhole. The evaluated action 



for each solution differs only in the sign taken for the square root in \\J a 4 — k 2 . The 
solution which is related to the Fourier transform of the wave function e~^ a2 cosh2 0-^-) 
corresponds to taking less than half the wormhole, and the positive square root. This 
means that the wave function goes like e - " 2 / 2 as a — > oo. The (a, 4>) representation has 
the advantage that the path integral appears to pick out the semi-classical wave function 
uniquely whereas the (a, k) representation does not. This is similar to the situation en- 
countered by Hartle and Hawking in defining the semi-classical approximation to the wave 
function of the universe in pure Einstein gravity with a cosmological constant [18]. In that 
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case, there are two euclidean four-geometries satisfying the field equations with a three- 
sphere boundary of given radius a a and no other boundary, but if instead the conjugate 
momentum to a is specified on the boundary, then there is only one. 

The analysis above shows a connection between the Fourier transformed wave function 
V> and the imaginary scalar field wormholes. We see from it that it is natural to take 
the boundary condition that the wave functions fall off like e - " 2 / 2 at large a. These 
wave functions correspond to saddle-point evaluations of path integrals where the four- 
geometries are taken to be asymptotically flat. We obtained a representation of the states 
\4>-) in terms of either (a, 0) or (a, k) . The (a, k) representation of can be related 
to imaginary scalar field wormhole solutions of the euclidean field equations. The (a, 0) 
representation suggested the correct boundary conditions for wormhole wave functions at 
large a . The next step is to ask whether it is possible to use the wave functions we have 
found to evaluate the vertex operators, and thus the effective action, corresponding to 
the euclidean wormhole solutions. We proceed by analogy with the conformally coupled 
scalar field [3] and the case of massless fermions, photons or gravitons [6-8]. Namely, we 
assume that there is a Hilbert space spanned by wormhole quantum states, \4>~) ■ The 
states \k) = J d(p- e ~ lfc< ^- also span this space. Each state has a vertex operator 

associated to it, which reproduces the effect of the state on low-energy n -point functions 
in the asymptotic region. One can try to calculate the vertex operators corresponding to 
the basis , where one fixes (j) — > 4>- at infinity. But the vertex operators turn out to 
be especially simple if the basis \k) is used. In the (a, 4>) representation the wave functions 
for these states are 'plane wave' products, of the form 

MaA)=K_ ¥k (a 2 /2)e-^ . (2.11) 

These wave functions satisfy the Wheeler-deWitt equation and are exponentially damped 
at large a . This boundary condition is suggested by the previous discussion of the saddle- 
point approximation to a path integral representation of the relevant wave functions. We 
are, however, being more general than [4] in our choice of boundary conditions in that 
we do not require any particular regularity condition at small a. Our view is that the 
Wheeler-deWitt equation is only an effective theory of quantum gravity, which ceases to 
be strictly valid at small a. (The semi-classical approximation breaks down there.) Thus 
any boundary conditions at small a should come from a more fundamental theory of 
quantum gravity, and not be imposed without justification. 

The effect of these wormhole states on the field theory n -point functions is given by 
the matrix element 

(k\(p(xi) ■ ■ -(j)(x n )\0) , (2.12) 
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as in [3] . This is essentially the dilute wormhole approximation, in that we are considering 
each wormhole independently and assuming that Green functions with points in the two 
asymptotically flat regions factorize. (We suspect that wave functions would have to be 
replaced by density matrices in a description which goes beyond the dilute wormhole 
approximation [19].) The vacuum state |0) is given by |0_ = 0) . The matrix element has 
a path integral representation as 

Jdaoij(ao) J dcPoMaoAo) J[dg][d<P]e- I ^^ ( j ) (xi)---<P(xn) , (2.13) 

where the inner path integral is taken over asymptotically flat four-geometries with an 
inner boundary on which a = cio and = 0o , and with <fi — > at infinity. We integrate 
in the range 0<ao<oo, — oo < </>o < oo with some unknown measure fJ,(ao) . The 
saddle-point approximation, applied to the inner integral, produces the result 



Jd^x Jda A(a ) J dfo ^(a , <fo) e -^ a o co8h2 ^ Qajj sinh20 o ) n^'-^o)" 2 , 

(2.14) 

for the matrix element, where the new co-ordinates x are related to the old t co-ordinate 
by 

(x-x ) 2 = -2t , (2.15) 

and we use (2.4) -(2.6). The co-ordinates x M become the usual Cartesian co-ordinates in 
the asymptotically flat region, and xq can be interpreted as the position of the wormhole 
end in asymptotically flat space. To obtain (2.14), we have also expanded the logarithm 
which occurs in the saddle-point solution (2.4) for 0, in the asymptotic region tj — > — oo. 
The prefactor A(a ) comes from the determinant of the fluctuations about the saddle- 
point, and from fi(ao) . Its exact form need not concern us, and is irrelevant in the leading 
semi-classical approximation. The scalar field action being purely quadratic, we know 
that A is just a function of a . The variables x are the zero-mode co-ordinates. They 
must be integrated over at the end of the calculation. In the simple model in which (p is 
homogeneous on the three-sphere spatial sections there is no additional integration over 
angular zero modes although in a more general model these integrals would also be present. 

Now the integral over </> may also be approximated using the saddle-point method. 
(This is an ordinary integral over <p Q in the spatially homogeneous model.) The saddle- 
point values for (p are the roots of 

ik = a 2 sinh 2</>o sp , (2-16) 
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so that in this approximation each of the n factors of ^aosinh2(/>o are replaced by their 
saddle-point values, \ik. (There is also another prefactor, which depends on k and ao, 
but not on n.) The ao integral gives a factor which will depend on k but not on n. We 
shall therefore ignore it, as it can be absorbed in the normalization of the vertex operator. 
The matrix element in the semi-classical limit (and for \xj — xq\ — > oo) therefore assumes 
the form 

J d 4 x a(k)(^ik^ flixj-xo)- 2 . (2.17) 

This can be identified with a flat space correlation function (V k (f>(xi) ■ ■ ■ (f)(x n )) , pro- 
vided the operator V k is taken to be 



yk 
v n 



! d 4 x Q ^ : (ik<P(x Q )r: . (2.18) 
J n\ 



Thus we obtain that the semi-classical effect of the wave function ip k (a 1 (p) on n-point 
functions in the asymptotic region is the same as the insertion of the operator V% in the 
flat space n-point functions, for some function a(k) . Crucial to this identification is the 
fact that the flat space (f> propagator, ((f) (x) (f)(0)) = \x~ 2 , falls off with the same power 
of proper distance in the asymptotic region as the saddle-point solution for <f)(x) , given by 
(2.4) with </>_ = 0. The factor of n! in (2.18) is present because each of the n factors in 
the vertex operator can be contracted with any of the n fields in the asymptotic region, 
and there are n! ways of doing this. The vertex operator we have found depends on n, 
but we require a vertex operator which reproduces the matrix element for any n. It must 
take the form of a sum 

V k = Vn = / d4x o a(k) : e lk ^ Xo) : . (2.19) 
o 

The effect of this operator on n-point functions is the same as (2.18) , to lowest order in 
h, which can be reinstated by multiplying each contraction of 0's by h, and dividing ik<f) 
in the vertex operator by h . The correct result to lowest order in % is all we can expect to 
obtain, since we have used the saddle-point approximation throughout. At this point we 
would like to remark that this analysis has just taken into account a single wormhole end. 
The inclusion of an arbitrary number of wormholes will lead to an effective integral over 
the coupling constants a(k) , for each value of k. This is similar to a functional integral 
over the field a(k) on superspace. a becomes a kind of quantized field on superspace 
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[9]. Whether there is some mechanism such as 'the Big Fix' [2] to fix the values of a(k) 
uniquely, remains unclear. 

The vertex operator (2.19) is exactly what was found in [16], which uses the imag- 
inary scalar field wormhole solution. However, it is important that nowhere here is it 
necessary to assume the existence of a complete wormhole solution interpolating between 
two asymptotically flat regions. We believe that, being not so reliant on exact wormhole 
solutions to the classical field equations, this derivation of the effective interaction will 
generalize more readily to more realistic models, such as those in which the restriction of 
homogeneity is removed. We shall now go on to discuss these more general models. 

3 Spatially inhomogeneous perturbations 

It has been suggested by Grinstein et al [5] that the vertex operator for massless scalar 
field wormholes given by (2.19) is correct, to all orders in l p i, the Planck length, and to 
lowest order in % . The calculations of [5] are based on the Green functions of test fields in 
an O (4) -symmetric euclidean wormhole background. We shall see in this section that to 
obtain derivative effective interactions (which will be higher order in l p i), it is necessary 
to consider inhomogeneities on the three-sphere sections. There will not be a complete 
wormhole solution connecting two asymptotically flat regions, but in the wave function 
approach this does not matter. What is important is that there is a family of solutions 
to the Wheeler-deWitt equation with appropriate boundary conditions. We shall start by 
investigating the wave functions in the midi-superspace approximation. Then the general 
form of the wave functions and the asymptotic behaviour of the modes will suggest the 
form that the effective interactions take. 

A The Wave Functions 

We shall assume a metric of the form 

ds 2 = a 2 [(N 2 + NiN i )dr 2 + 2N t dx l dr + h i3 dx l dx j ] , (3.1) 

where the three-metric hij is 

h ij = e 2a ^(il ij + e ij ) . (3.2) 
Here Vtij is the round metric on the unit three-sphere, and we expand e^- in harmonics, 

Cij = ^2 W^ a riirn-^Q 7 lrrS^ij + ^^nlm (PijYlm + ^ c nim (<% )™ m + 2d n £ m (Gij)^ m \ , (3.3) 
n£m 



10 



with the sum starting at n = 2, and the notation of [10]. The n = 2 modes are rather 
special: the n = 2 traceless tensor harmonics all vanish identically and so the n = 2 term 
in the sum is just the ^a2Q 2 Slij part. (We shall henceforth drop the labels £m.) The 
lapse, shift and scalar field are similarly expanded: 



N = Ni 



! + 6- 1/2 E^ r 



N i = e a J2i 6 ~ 1/2 ^ + (St) n ] 



(3.4) 



$ = a" 1 



V , rl 



The action is expanded to all orders in a, 0, No and to second order in the perturbation 
quantities q n G (a n , 6 n , c n , <i n , /n) and r n G (g n , k n , j n ) . We refer the reader to [10] for 
the details of the action, constraints and three-sphere harmonics. Conjugate momenta 
can be defined in the usual manner, and the hamiltonian expressed as a function of co- 
ordinates and momenta. (We denote the euclidean momenta conjugate to q by ir q .) The 
hamiltonian takes the form 



H = N, 



' H n A- i V 



ir\) 



(3.5) 



with the subscripts denoting the order of each part of the hamiltonian in the perturbation 
quantities, and whether each part arises from varying the lapse or shift. 
We look for wave functions which can be expressed in the form 



(3.6) 



where ^o(ct, <j>) is one °f the euclidean WKB solutions of the mini-superspace background 
Wheeler-deWitt equation discussed in Sec. 2 (for instance one of the ipk(e a ,(fi) defined 
in (2.11)). The full Wheeler-deWitt equation (H\ + En ^2)^ = reduces to a set of 
euclidean Schrodinger equations for ifj n , along the euclidean trajectories defined by the 
background wave function. The partial wave functions tfj n must also satisfy the linearized 
hamiltonian and momentum constraints. The euclidean trajectories which are defined by 
the background wave functions ipk are the solutions to the background classical equations 
of motion (2.4) , with Q = ik. We shall consider this choice of background wave function 
from now on. The vertex operator corresponding to this choice will contain a factor of 
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e tk<j} ^ f rom 1 p k ? anc [ a factor from each of the ip n , which will depend on which solution to 
the euclidean Schrodinger equation is chosen for each n mode. This framework, with a 
WKB background plus quantized perturbations satisfying a Schrodinger-type equation, is 
the 1 M -expansion', which is discussed further in [20]. We shall later discuss qualitatively 
the type of wave functions one might expect from the modes with n > 2, but first we 
consider the n = 2 modes. (We drop the mode label 2 on the perturbation quantities for 
ease of notation.) 

The partial wave function ifj 2 depends just on a, /c, a, /. The n = 2 linear constraint 
equations to leading order in the M -expansion are 

S H_iip2 = V 3a (-7r a + an a + 3/tt^) ^ 2 = 

: (3.7) 

#|iV>2 = -e" 3a {a(nl + 3^ - 3e 4a ) - 2(n <p n f - 7r a 7r a )) ^ 2 = . 

The carets on 7r ,7r/ denote the operator form of the euclidean momenta conjugate to 
a, / (represented as —d/da,—d/df), while the quantities 7r a and 7r^ are the euclidean 
momenta of a, 4> in the background solution. These are given by ix a = \/e 4a — k 2 and 
7r^ = ik. We can use the first of equations (3.7) to substitute for n a ip2 into the second 
equation, which we can then solve on the surface a = 0. This can then be used as an 
initial condition for the first equation, to obtain ^(a, k, a, /) at non-zero values of a. We 
see that ip2 has the form of a Gaussian at a = , 

M^k,0J) = c 2 (a,k)e- L ^ a ^ f2 , (3.8) 

although it is not clear what i\)<i represents since both a and / are gauge quantities. The 
main point is that the linearized constraint equations have picked out one solution to the 
euclidean Schrodinger equation for these modes, which appears to be in a kind of 'ground 
state'. 

For n > 2 the situation is completely different from the above. In this case there 
are three linearized constraint equations for five perturbation quantities, so there are two 
physical degrees of freedom remaining; one scalar (say s n , originating from a n , b n , f n ) and 
one tensor (d n ). We can perform a reduction to the physical modes before quantization 
(see [21,22]). Then each of the n > 2 partial wave functions ip n depends on a, k, s n , d n 
with the s n and d n parts being separable. Each part obeys a euclidean Schrodinger 
equation of the form 

N P h H^ n = -—i; n , (3.9) 
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where the operator ph H, 



is a diagonalizable homogeneous quadratic in q n and jt qn , and 



Qn £ (s n ,d n ). The problem thus reduces to the euclidean Schrodinger equation for a 
harmonic oscillator with time-dependent frequency. (For the tensor modes d n in pure 
gravity this equation reduces to the simple harmonic oscillator, as shown in [8].) The 
method of Salusti and Zirilli [23] (also discussed in [24]) can be used to find a complete set 
of orthonormal solutions. We describe how this can be done, and the inner product defined 
in the euclidean sector, in the Appendix. The states can be obtained from a Gaussian 
ground state by the application of suitable 'raising operators', in a similar manner to the 
simple harmonic oscillator. The action of m raising operators on the ground state yields 
a state which is interpreted as describing a closed universe containing m particles in that 
mode. It takes the form 



where H m is an mth order Hermite polynomial. The solutions to the euclidean Schro- 
dinger equation along the classical trajectories of fixed k can then be 'lifted' to full solutions 
of the Wheeler-deWitt equation, using the method described in [20] . 

B The Classical Solutions and the Vertex Operators 

There are two ingredients in the path integral formula for the matrix element (2.12) . 
One is the wave function, which describes the quantum state of the wormhole. We have 
described above the general form that such wave functions take. The other is the path 
integral over asymptotically flat four-geometries with an inner boundary on which the 
three-geometry and matter field are specified. For the semi-classical evaluation of this 
path integral we need the classical euclidean solutions for the perturbation modes. We 
shall now discuss these. 

The classical euclidean solutions including the n = 2 modes can be obtained from the 
background solution (2.4) by slicing the background in a different way. The background 
solution can be expressed in the form (oo > r > r > y/\Q\/2) 



l(>nm(T, q n ) (T)H m (p n (r)q n )e 



(3.10) 



ds 2 = a 2 n(r 2 )(dr 2 + r 2 dn 2 3 ) 
<P = <P(r 2 ) , 



(3.11) 



where 




-2 




(3.12) 
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We now change co-ordinates from r, x to r',x', with r cos x = r' cos x' — e and r sin x = 
r'sinx' (so that r 2 = r' 2 — 2er'cosx' + e 2 ), for a constant parameter e. This change is 
an isometry of the flat metric dr 2 + r 2 dfl 2 . However O and (p develop non-spherically 
symmetric parts in the new slicing of R x S 3 in three-spheres of constant r' . (We also 
change the boundary to one at constant r' , instead of at constant r.) We see that 

0(r 2 ) = fi(r' 2 ) - eQ 2 r'- 5 cos X ' + 0(e 2 ) 

(3.13) 

(j)(r 2 ) = ( j)(r' 2 ) + eQr'- 3 n- 1 (r' 2 ) cos X ' + 0(e 2 ) . 

as e — > . We thereby obtain an exact solution to the linearized equations for the n = 2 
modes, since a/2/7t cos % is an n = 2 scalar harmonic. We shall now remove the primes 
on r',x' • (Notice that in this form the solution is obtained in the particular gauge k = 
0, g = a, N = e a r~ 1 .) 

The n = 2 scalar field perturbation falls off like r~ 3 cosx in the asymptotic region 
(r — > oo), while the background part of the scalar field falls off like r~ 2 . It is essentially 
the fact that (p falls off with proper distance like the two-point function in flat space which 
led to the vertex operator being a function of (p (no derivatives of <p) for the spherically 
symmetric model of Sec. 2. We see from (3.13) that the n = 2 mode falls off like a single 
derivative of the two-point function c^d^r~ 2 = — 2r -3 cosx, where c M is a constant vector 
in flat space. Thus the vertex operator coming from these modes would be expected to be 
a function of <9 M . By considering the asymptotic form of the linearized equations for the 
general n mode, we see that the scalar part falls off like r~ n ~ 1 Q n , for a scalar harmonic 
Q n . Thus the contribution from the nth mode will produce a vertex operator which is 
made up of n — 1 derivatives of <p. There is an integration over the 0(4) rotations at the 
end of the calculation (a zero mode integration, similar to the integral over xq in Sec. 2). 
This means that the total vertex operator will be Lorentz covariant. In general, a wave 
function in the scalar sector of the form (3.10) would be expected to lead to a vertex 
operator (V n ) m which is the mth power of a basic operator V n . The operator V n will be 
made up of n — 1 derivatives of (p . However, we have found for the n = 2 mode, that the 
wave function satisfying the linearized hamiltonian and momentum constraints is forced 
to have Gaussian form. For such wave functions the vertex operator will be restricted to 
have m = 0. Thus we have appeared to rule out the operators {defy m , for m > 0, but 
not corresponding operators involving higher derivatives of (p. This appears to contradict 
work in [11], which considers the conformally coupled massless scalar field n = 2 mode. 
In [11] it appears that the vertex operator (dep) 2 is present, corresponding to the lowest 
allowed excited solutions to the Wheeler-deWitt equation. We expect however, that the 
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conformally and minimally coupled massless scalar fields should have the same number of 
degrees of freedom for each level n , and that our result on the absence of first derivative 
terms in the vertex operator would apply also in the conformally coupled case. 

Our argument fixing the form of the n = 2 wave function to be in its ground state 
relied on the fact the n = 2 modes are pure gauge. It does not apply to the modes n > 2 . 
For the n > 2 modes excited states will exist which satisfy the Wheeler-deWitt equation 
and linearized constraints. It is expected that the part of the vertex operator (in the 
scalar part of the theory) coming from the mth excited state will consist of m powers 
of operators involving n — 1 derivatives of 4>. This will yield corrections to the vertex 
operator e lk ^ in the massless scalar field model which cannot be obtained using the Green 
functions of fields in the spherically symmetric wormhole background [5]. Thus the issue 
of whether the semi-classical vertex operator has really been found to be e lk ^ to all orders 
in the wormhole scale, is still a matter of debate. 



4 Conclusions 



We have discussed the form of the vertex operator for wormholes in the massless mini- 
mally coupled scalar field model. We used the wave function approach of Hawking [3], and 
described the evaluation of vertex operators corresponding to a basis of solutions to the 
hamiltonian and momentum constraints. The O (4) -symmetric mini-superspace truncation 
yielded a correspondence between a particular basis of wave functions which separated in 
the (a, 4>) mini-superspace variables, and effective semi-classical vertex operators of the 
form J d 4 x a(k): e lk ^ x °^>: . These vertex operators had previously been obtained by con- 
sidering the field theory in the background of an imaginary scalar field euclidean wormhole 
solution [16]. Using that method it appears that there are no further corrections to the 
vertex operator of higher order in l v \ (or the wormhole characteristic scale). However, 
using instead the wave function approach it appears that such corrections exist, and would 
come from non-spherically symmetric perturbation modes on the spherically symmetric 
background. The n = 2 modes, which would yield a factor in the vertex operator like 
(d(/)) m , appear to be pure gauge, and so the linearized hamiltonian and momentum con- 
straints restrict this part of the wave function to be in its ground state. Thus these terms 
appear not to be present. This appears to contradict earlier work on the conformally cou- 
pled scalar field model [11]. However the n > 2 modes can be excited, and would appear 
to lead to higher derivatives in the effective interaction. This contrasts with the results ob- 
tained in [5] using the field theory in the imaginary scalar field wormhole background, and 
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suggests that the two different approaches to wormholes correspond to physically different 
situations. 
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Appendix 

We show here how to construct the basis of solutions to the euclidean Schrodinger 
equation which are used in Sec. 3. We also construct an inner product which is particularly 
natural for wormholes, and with respect to which the wave functions are orthonormal. 

Consider the case where the euclidean WKB background wave functions are the 
ip k (e°' 1 (p) defined in (2.11). Each wave function ip k defines a euclidean Hamilton-Jacobi 
function 

I k (a,(f>) = i (Ve 4Q - £; 2 + £;arcsin(£;e- 2a )) + ik(p . (A.l) 



We shall denote the partial derivatives by I a = Ve 4a — k 2 and 1$ = ik . There are 
associated euclidean trajectories, which are the solutions to the equation of motion (2.4) , 
with Q = ik . The classical euclidean momenta are given by n a = J Q , = 1^. Since the 
background wave functions have WKB form, in the region e 4a ^> k 2 ^> 1 , we can reduce 
the quantization of the n > 2 perturbations to the study of the Schrodinger equation for 
one tensor and one scalar mode, in the background solution [22]. These are the physical 
degrees of freedom, once all linearized constraints have been solved. In our case, the 
background classical geometries are euclidean, and one obtains a euclidean Schrodinger, 
or heat, equation. The equation is generally of the form 

1 d 

H n ifj k (T,q n ) = - — —ifj k (T,q n ) , (A.2) 

where the 'time parameters' r along the trajectories have been introduced via <fre 3ct = NqI^ 
and — cte 3a = NqI^ . This can be conveniently solved in the gauge Nq = e~ a by r = ■ 
Leaving out the suffices n, and choosing a convenient operator ordering, the 'hamil- 
tonian' H is 

H = l -e~^ (a% + l -B{qd q + d q q) + Cq 2 ^ , (A.3) 
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where, for the tensor modes (q = d), A = — 1, B = 0, C = e 4a (n 2 — 1) , and for the scalar 
modes (q = s = f + ^(a + b)), 



A = - 1 + 



B = Q( r ^\ I i (A.4) 




n 2 - 4 J / u 

c = <" 2 - d ('« - 

The same procedure as used in [22] has been used to obtain these equations, except for 
the difference that this analysis is in euclidean time. 
By writing 

ip k (r,q) =e~^ q2 Xk(r,q) , (A.5) 
the equation for Xk takes the form of (A. 2) with A —> A, B — > and 

C~C> = C . (A.6) 



AA 2N dr \A 
Solutions exist for Xk which are of Gaussian form 



for any u satisfying the linear equation 

1 d fiie 3a 



N Q dr \N A 



+ C'e- 3a u = . (A.8) 



The equation (A.8) has a number of important properties. Firstly, it is real and so it 
possesses two linearly independent real solutions. This is not the case for its lorentzian 
counterpart. Secondly, in the asymptotic region, t ^ — oo, A ^ — 1 and C ~ e 4a (n 2 — 1) , 
so (A.8) reduces to the equation for the tensor modes. One of the solutions to (A.8) 
tends to zero while the other tends to infinity, in this region. We can use this to define 
a preferred vacuum wave function, Xko(T,q) which is obtained by choosing the solution 
to (A.8), u(r) = u (t) , which tends to zero in the asymptotic region. (In the gauge 
iVo = 1 we have uq(t) ~ r~ n_1 .) Thirdly, we note that (A.8) is unchanged with respect 
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to t — > —t . (If we choose, for instance, the gauge N = e~ a , with the choice of time 
parameter r = — |/ Q , this is made somewhat more explicit.) The origin of this symmetry 
is the discrete isometry of the wormhole geometry, which interchanges the two asymptotic 
regions. This symmetry means that if u(t) is a solution to (A. 8) , then so is u(t) = u(—t) . 
The 'vacuum' wave function, xw , is annihilated by the operator 



ue 3a 



a(u) = ud q - jj-jq , (A.9) 

with u = uq . The operator a(uo) , when applied to xw , yields another solution to the 
euclidean Schrodinger equation, (non- vanishing, because uq(t) and uo(r) are linearly 
independent). We call this solution Xki- This process can be continued indefinitely to 
produce Xkm- The operator a(uo) is the 'generalized raising operator'. We can then 
construct ipkm from Xkm using (A. 5) . 

To construct the appropriate inner product, in which these wave functions are all 
orthogonal, we first use uq(t) to define a 'dual vacuum' 



(A.10) 



~ / x —1/2 / u e 3a 2 \ 
Xko(r, q )=u expj^-^-^, j 

ipko(r,q) = e + ^ q Xko(r,q) , 
which satisfies the 'time reversed' euclidean Schrodinger equation 

\e~ Za (A% - \B(qd q + d q q) + Cq 2 ^ &(t, q) = + l-£-j, k ( T , q) . (A.ll) 

(Under time-reversal we note that B — > —B.) The operator 

a{u) = -ud q - j^-^q , (A.12) 

annihilates x^o , if u = uq , and if u = uq it can be applied m times to Xko , to produce 
Xkm and hence the 'dual mth excited state', ifikm- 
With these definitions, the inner product 

/+oo 
rf#(r, ? )V(r,Q) (A.13) 
-oo 

is independent of r, and with suitable normalization constants the wave functions ipkm 
are ortho normal with respect to this inner product. 
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The wave functions V>fc m (r, Qn) have the form quoted in the text, 

ipkm{r,q n ) =-fknm(T)H m ((3 kn (T)q n )e- akn{T)q * , (A.14) 

where H m is an mth order Hermite polynomial. They are interpreted in terms of quantum 
states corresponding to a closed universe containing m particles in the nth harmonic mode. 
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